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1. Introduction 

The Baxter's three-color problem on a regular hexagonal lattice is one of the 
classical examples of solvable lattice models. The free energy of this model is equal to 
the number of coloring the links of the lattice with three different colors, A, B and C, 
so that no links that meet at a vertex carry the same color. Alternatively one can speak 
of the number of three-coloring of the links of the regular triangulated lattice in a way 
that the three sides of every triangle have different colors. The problem is equivalent to 
a special version of the 0(2) model, whose partition function is given by a gas of fully 
packed loops on the hexagonal lattice, having two different flavors. Baxter also showed 
that this model solves the problem of coloring the faces of the hexagonal lattice with four 
different colors, so that the adjacent faces have different colors. More recently yet another 
geometrical interpretation of this model has been found, namely as the problem of counting 
the different foldings of a regular triangular lattice ||. 

The three-coloring problem can be formulated also for a random 3-coordinated planar 
graph (or a random triangulation, in the dual language) 0. Here one has the freedom 
to chose different ensembles of planar graphs. The generating function for all possible 
three-colored planar graphs is given by the three-matrix integral 0] 

Z N {j3) = J dAdBdC exp!^-NTr(^(A 2 + B 2 + C 2 )-g[ABC + BAC}^ (1.1) 

where A, B and C are hermitian NxN matrix variables. 
The perturbative expansion of the free energy 
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has the following combinatorial meaning: Nh(2rn) is the number of the connected three- 
colored triangulations of a surface of genus h, containing 2m triangles. (The configurations 
with nontrivial symmetry group are taken with the corresponding symmetry factors.) 

The model has still the interpretation as a fully packed 0(2) loop gas model, with 
the restriction that only loops of even length are allowed. It can be also interpreted as the 
four-coloring problem of the faces of a random 3-coordinated graph. The interpretation 
in terms of foldings does not generalize to this ensemble: it holds only for planar graphs 
whose faces have even number of sides. 

An iterative procedure for evaluating the coefficients in the expansion (|1.2|) has been 
proposed by B. Eynard and C. Kristjansen ||. They considered, more generally, the 0(2n) 
loop gas model with only even loops, which reduces to the three-color problem when n = 1. 
This 0(2n) model is actually identical to fully packed 0{n) loop gas model defined on a 
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4-coordinated planar grapbM It is therefore expected that the matrix model ( |1 . 1|) is in the 
class of universality of the dense 0(1) model, which means that A/) l (2m) have the same 
large-m asymptotics as the number of the non-colored planar graphs: 

N h {2m) ~ g*~ 2m m^ h -V-\ (1.3) 

where g* is the critical value of the coupling g. The iterative procedure was carried out up 
to the order 12, with the result 

lim F(g) = 2g 2 /2+lAg A /A + 138# 6 /6+1608# 8 /8+20736# 10 /10+286452# 12 /12+... (1.4) 

Soon after the work of B. Eynard and C. Kristjansen, the exact solution of the matrix 
model ( |L2^ ) has been found for purely imaginary coupling g || . An explicit expression for 
the free energy of the model was obtained in terms of elliptic functions, and it was checked 
that the coefficients in the expansion (|1.4j ) are correctly reproduced However, it was 
not possible to to carry out the analytic continuation of this solution to real values of g. 

In this paper we present the direct solution of the matrix model ( |1.2|) for real coupling 
g, using a method similar to the one applied in ||. We find the expected critical behavior 
of a c = matter coupled to 2D gravity and evaluate the critical coupling g*. Finally 
we we explain why the cases of real and imaginary coupling g are not related by analytic 
continuation. 



2. Saddle point equations for the matrix integral 

After integrating over the B and C matrices, the partition function reduces to the 
following integral over the eigenvalues Ai, . . . , Ajv of the matrix A: 



The saddle-point spectral density 
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(2.2) 



is an even function supported by a symmetric interval [—A, A], where A is a function of 
the coupling g, such that A| g= o = 2. At g = it is given by the Wigner's semi-circle law: 
p(X)\ g =o = t~ V4 — A 2 . It is completely determined by the saddle-point equation 



A = -IdX' p(X') 

J-A 
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A -A' X + X' + g- 1 X + X'-g- 1 



(2.3) 



1 More precisely, each 4-vertex is visited either by one loop going straight, or by two loops 
turning at right angle. 
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and the normalization condition J dX p(A) = 1. Introducing the analytic function 

W(z) = f A (2.4) 

we write the saddle point equation as a functional condition for W(z): 

W(X + iO) +W(X-iO) + W(-l/g - X) + W(l/g - X) = A, Ae[-A,A]. (2.5) 

The density p(A) is an even function of A, which means that (|2.4|) is antisymmetric: 

W(z) = -W(-z). (2.6) 

The saddle point equation can be therefore rewritten as 

W{X + iO) + W(X-iO)-W(X + l/g)-W(X-l/g) = X, Ae[-A,A]. (2.7) 

The normalization condition for the density fixes the first coefficient in the expansion of 
W(z) at infinity: 

W(z) = -+ 1 +0(z 5 ). (2.8) 

z z^ 

As mentioned in ||, knowing the the quantity (j^-A 2 ) is enough to solve the three- 
coloring problem, because it is related to the free energy by 

w A2 ) = 1+ %™ < 2 ' 9 > 

3. Solution of the saddle point equation 



The four-term functional equation (277) with the boundary condition W(z) ~ 1/z at 



infinity is sufficient to determine W(z). Introducing the function 

C(z) = 2g- 1 [W(z + 1/20) - W(z - l/2g)} + z 2 , (3.1) 
we write ( |2.7| ) two-term difference equation 

C(A + l/2g ± iO) = C(A - l/2g T iO), A G [-A, A]. (3.2) 
The function ((z) is real and symmetric: 

c(z) = a- z )=ai), (3.3) 
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possesses two cuts along the intervals [—A — A — and [—A + A + i] on the 
real axis, and expands at z — > +00 as 



C(*) = * 2 - 



(^+6^ 2 ) 



g 2 z 2 



9 z* 



+ 0(z~ 6 ). 



(3.4) 



From the symmetry of ((z) it follows that it is real when z is real and outside the 
cuts, and also when z is purely imaginary. Therefore the map z — > C transforms the 
quadrant {Imz > 0, Rez > 0} into the to the upper semi-plane {Re£ > 0}, cut along 
an arc connecting the points c and b, where b is some complex number (Fig. 1). At the 
endpoint c, the arc should be perpendicular to the real axis. We denote the images of the 
special points of the map 



by 



1 A 1 1 A 

zi=0, z 2 = - A, z 3 = — , z 4 = — + A, 

2g 2g 2g 



a = C(zi), b = C(z 3 ), c = C(z 2 ) = C(z4). 



(3.5) 



(3.6) 





l/2g- A l/2g 

Fig.l: The map z <-> C- 



l/2g+A 



It is easy to see that the inverse function z(C) can be obtained as the integral^ 

1 r c dt (t - c) 



2Ja y/(t-a)(t-b)(t-b)' 
The asymptotics ( ^.4| ) is satisfied if 

a + b + b - 2c = 0, 2c 2 - a 2 - b 2 - b 2 = 6/g 2 . 



(3.7) 



(3.f 



Here we are following the argument by J. Hoppe in [^]. 



The integral ( |3.7|) is a standard elliptic integral. Denoting 



cnu = cn(u, k) = - — b - a = Ae ie , = cos(0/2). (3.9) 



A + C-a 

we have ( @, 239.00, 239.07 and 341.53) 

z a — c — A 



y/A 
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du 



o 1 + cnu 



where 



Cu - H ' {U ' k) + dn{u ' k) 
H(u, k) sn(w, k) 

A+a-c E 

O — 



(3.10) 



2A K- < 3 - U » 
The elliptic parameters u±, . . . , U4 of the special points of the map z — > £ are (see Fig. 2) 



Mi =0, -u 2 G [0,-fT], u 3 = K + iK',U4 = 2u 3 - u 2 , 
and the point z —* 00 (£ — > 00) corresponds to Uqo = 2zi^ ; . 



(3.12) 
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2K 



Fig. 2: The domain of the u- variable. 



From ( |3.5| ) we determine the constants A and C 
z 3 = z(K + iK') = VA 

which gives 



7T 



-C.h' + / ( — + CK' 



C 



n 



Va=^k'. 



1 

2? 



2KK' 7T 

Then, using the Legendre's relation EK' + E'K — KK' = 7r/2, we write ( 3.11 ) as 

c - a _ 2E' - K' 
A ~ ~K' ' 
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(3.13) 



(3.14) 



(3.15) 



The final formula is 



. z(u) tv H'(u,k) dn(u,k) 

Z[U)= ~7t = ZKK' U + ~H(^k) ~ M^k) 

~ _ ((u) - a _ 1 - cn(u, k) 

;W A l + cn(u,fc)" 



Note the very useful relation 



dz ~, . 
du 



2E' - K' 



(3.17) 



Now it remains to determine the elliptic modulus as a function of g. It is fixed by the 
Jj-term of the expansion 



2 1 ^ + 6<^4 2 >1 



A g 2 A 2 z 2 



g 2 A 3 



+ 



Comparing the Laurent expansions of ((u) and z(u) at the point Mqo 
expansion of z(u) is most easily obtained using the relation (|3.18|)), we find 
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7T 

~K' 



IE' -K' 1 - 2k 2 
2K' + 3 



1 - 16k 2 k' 2 
36 



(3.18) 
UK' (the 

(3.19) 



The expansion ( p. . 2|) of the free energy can be obtained using (|2.9| ) from the ^-term of 
the expansion ( |3.18| ). The Gaussian limit g — > corresponds to k' = 4g + g 2 + ...—> 0. 
Therefore, it is convenient first to express the solution (|3.16|) in terms of the dual modulus. 
Using the duality relations 

H(u,k) = -i^KjK'e-^l^' H(iu,k') 
1 



cn(u, k) 



ds(w, k) = ids(iu, k') 



we write 



z(u) = i 

c» = 



cn(m, k') 

H'(iu,k') dn(iu,k') 



H(iu, k') sn(m, k') 
cn(zw, /c') — 1 
cnfzw, fc') + 1 



(3.20) 



The large- AT limit of the free energy T is obtained from the asymptotics (|3.18| ). This is 
purely technical exercise and we will not do it here. Instead, we will obtain the asymptotical 
behavior of the free energy when the volume of the graph diverges. 
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4. Critical behavior 



The Taylor expansion of the free energy in g is convergent up to the critical point 
<7*, where the nearest singularity of J-(g) is located. Since the free energy is an analytic 
function of k 2 = 1 — k'J' , the critical point is determined by d/3/dk = 0. The function f3(k') 
has a minimum at ki = 0.826114 where 

2E(K)=K(K). (4.1) 

In the vicinity of the critical point (k' = + 8k') 

$r 1 - l ^W Lsk ' 2+0{Sk ' i) ' (42) 



where 



g> = KHK1 = (*")«-* . (4.3) 
7i 2 \/6 7r 2 -\/6 



By ( |3.15|) this is exactly the point where the the right end of the left cut of touches 
the left end of the right cut. At this point u\ = u\ = 0, c* = a* and A* = Finally, 
it is evident that the critical singularity of the free energy is F(g) — ^(g*) ~ (<?* — fif) 5 ^ 2 , 



which implies (1.3). 



5. Conclusion 

We have found that the coefficients of the free energy ( |1 . 2|) grow as 

A/- (2m)~^- 2 ™ m- 7 ' 2 . 



On the other hand it is known [|1(J that the number of 3-coordinated planar graph with 
2m vertices grows as (12\/3) n n~ 7 ^ 2 . Therefore the number of three-colorings per vertex 
of an infinite random 3-coordinated planar graph is equal to 

J/</: - 0.984318.... (5.1) 



For a large but finite planar graph, the number of three-colorings grows linearly with the 
number of vertices A = 2m, up to a term O(^). There is no logarithmic corrections. The 
fact that the number of three- colorings is slightly less than one is explained by the fact 
that not all 3-coordinated graphs are three-colorable. 



In terms of the 0(2) loop gas, eq. (|5.1|) gives the entropy per vertex of the gas 



fully packed loops on a random 3-coordinated graph, having two different flavors and even 



7 



length. This is to be compared with the entropy of the gas of fully packed loops with n 
flavors and no restriction on the length [11| 



\/g° w _ 2y/2(2 + n) 
\/\2y/% \/l2V3, 

which is equal to 1.24081 for n = 2. 

In order to interpret our result in term of the 0(1) loop gas on a 4-coordinated planar 
graph, we recall that the number of planar 4-coordinated graphs with m vertices grows as 
12 m m -7 / 2 ||10||. Therefore the entropy per vertex of the fully packed non-oriented loops 



on a 4-coordinated random graph is is given by 

= 1.62593.... 

12 

Finally, let us mention that the integral (|2.1|) is related to the Baxter's 6-vertex model 
on a random lattice, which has been recently solved exactly ]12] . In the parameterization 



a = b = 1, c = — 2 cos (3 
of the vertex weights, the integral (|2.1|) corresponds to the limit /3 — > 0, where (3 is assumed 



to be purely imaginary. The parameter 

a 2 + b 2 - c 2 



- cos 2(3 



2ab 

achieves the value —1 from the left, which means that we are approaching the boundary 
between the regimes III and IV from the side of the regime IV. This regime IV is char- 
acterized by a finite correlation length, which is confirmed by our solution. On the other 
hand, the matrix integral considered in || can be identified with the limit /? — ^ with 
[3 > 0, which means that the the boundary III/IV is achieved from the critical regime 
III characterized by an infinite correlation length. The fact that the real and imaginary 
values of the coupling g are associated with different regimes of the six-vertex model ex- 
plains why our solution cannot be obtained from the solution found in || by analytic 
continuation ig — > g. Even if the perturbative expansions in g 2 are identical, they differ 
by nonperturbative terms. 

Note added: After this manuscript was to the publisher, the author learned that Bertrand 
Eynard and Charlotte Kristjansen found independently the exact solution of the problem. 
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